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摘要 : 建立 了 以 无 数 个 无 限 深 势 阱 中 的 粒子 为 工 质 的 不 可 逆 量 子 斯 特 林 致 冷 循环 模型 。 势 阱 中 粒子 在 能 级 上 的 分 布 由 吉 布 
布 函数 决定 。 该 致 冷 循环 由 两 个 等 能 量 过 程 与 两 个 等 势 阱 壁 宽度 过 程 组 成 ， 具 有 高 低温 热源 间 的 热 漏 和 不 完全 回 热 因素 。 导 
了 性 能 系数 、 致 冷 率 、 烂 产 率 和 生态 学 目标 函数 性 能 解析 式 ， 制 冷 率 与 生态 学 函数 的 关系 曲线 是 回 原点 的 扭 叶 型 。 分 析 了 热 沁 
系数 与 不 完全 回 热 因子 对 性 能 系数 、 致 冷 率 和 生态 学 函数 的 影响 。 
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Study on the performance of quantum Stirling refrigeration cycle 
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Abstract: An irreversible quantum Stirling refrigeration model is established in which working substance is a single 
particle contained in a one-dimensional infinite potential well, and the system consists of countless replicas. The 
occupation probabilities in each eigenstates of the particle are expressed by the thermal equilibrium Gibbs 
distributions. The cycle is made of two iso-energy branches and two isochoric branches and with heat leakage 
between heat reservoirs and imperfect regeneration. Analytical expressions of COP, cooling rate, entropy generation 
rate and ecological function are derived. The curve between cooling rate and ecological function is a loop-shaped 
one. Effects of heat leakage and imperfect regeneration on COP, cooling rate and ecological function are analyzed. 
Keywords: finite time thermodynamics; one-dimensional infinite potential well; quantum refrigeration cycle; 
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0 引言 


在 微 纳米 尺度 09、 低 温 或 高 密度 条 件 下 ， 需 要 ”引入 了 旁 通 热 漏 。Feldmann 和 Kosloff?7 引 入 了 内 
考虑 热力 循环 工 质 的 量子 特性 。 应 用 有 限时 间 热 力 ”摩擦 系数 来 描述 量子 非 绝 热 现 象 。 吴 锋 等 C9 考 虑 热 
学 理论 和 方法 55, 可 以 建立 各 种 量子 循环 模型 0520， ” 漏 、 热 阻 以 及 有 限 速 率 传 热 对 不 可 道 量子 斯 特 林 玮 
研究 其 热力 性 质 与 性 能 参数 ， 优 化 其 运行 区 间 。 实 。 冷 机 性 能 的 影响 . 金 晓 昌 多 从 庄 经 济 学 角度 研究 了 
际 热机 中 ， 存 在 着 热 阻 、 热 漏 、 摩 擦 、 涡 流 、 惯 性 ” 量子 斯 特 林 制 冷 机 最 优 性 能 。 

效应 及 非 平衡 的 影响 ， 循 环 为 不 可 逆 循 环 吴 。 一 些 在 文献 [3,10,14,16,23,28,30] 的 基础 上 ， 本 文 将 建 
学 者 也 研究 了 不 可 首 因 素 对 量子 热力 循环 性 能 的 影 。 立 一 个 包含 高 低温 热源 间 的 热 漏 与 不 完全 回 热 的 不 
响 。 金 晓 昌 等 29 考 虑 高 、 低 温 热 源 间 的 弱 耦 联 作 用 可 逆 斯 特 林 制冷 机 循环 ， 其 工 质 为 无 数 个 囚禁 于 一 
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维 无 限 深 势 阱 中 的 粒子 。 导 出 循环 的 性 能 系数 、 致 
冷 率 、 粮 产 率 、 生 态 学 函数 ， 分 析 这 些 重要 的 性 能 
参数 与 热 漏 率 和 不 完全 回 热 因子 之 间 的 关系 。 


1 系统 的 量子 力学 描述 


本 文 研究 的 对 象 是 一 个 囚禁 于 宽度 为 工 的 一 维 
无 限 深 势 阱 中 的 粒子 。 其 定 态 波 函 数 满足 : 


2 2 Se (D) 
d“Y/dx +2mET/h’ =0 


根据 波 函 数 到 必须 满足 的 边界 条 件 WO)=0 和 
WL)=0， 可 求 出 波 函 数 为 : 


(xX) 2 (2) 
ns 


式 


9 是 归 一 化 的 本 征 函 数 : 
$2)=V2/ Lsinnzx/ D) G) 


系数 四 满 足 归 一 化 条 件 ; 


a 由 


式 中 pn=aw?， 为 概率 密度 ， 系 统 的 能 量 本 征 值 为 : 

E,=n A (2mL) 已 
式 中 m 和 分别 为 粒子 的 质量 与 量子 数 ， 哈 密 顿 量 
的 期 望 值 为 : 


= Dp,E, (6) 
n=l] 


经 典 循 环 中 ， 活 塞 在 气缸 中 做 往复 运动 ， 假 设 
势 阱 壁 像 经 典 循环 中 的 活塞 一 样 运动 , 当 势 阱 宽 i 度 
变化 时 ， 系 统 的 波 函 数 与 能 级 都 会 随 着 也 变化 。 定 
义 广 义 力 为 : 


07) 
六 二 —dW / dy, 


式 中 ，y 是 与 广义 力 7 对 应 的 广义 坐标 。 因 此 ， 作 
用 在 势 阱 壁 的 力 可 以 表示 成 : 


F=-dE/dL= 2 p,ln A h’ / (mL)] (8) 


n=l1 


2. 量子 斯 特 林 致 冷 循环 
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研究 内 禁 于 一 维 无 限 深 势 阱 中 的 粒子 ， 为 简单 
起 见 ， 只 考虑 两 个 能 级 。 循 环 工 质 为 无 数 个 这 样 的 


粒子 ， 每 一 个 都 被 办 禁 在 一 个 一 维 无 限 深 势 阱 中 ， 
粒子 在 能 级 上 的 占有 几率 由 吉 布 斯 分 布 决定 。 循 环 


由 两 个 等 能 量 过 程 与 两 个 等 势 阱 壁 的 回 热 过 程 组 
， 分 别 与 经 典 循环 中 的 等 温 过 程 与 等 容 过 程 对 应 。 
考虑 不 完全 回 热 与 高 低温 热源 之 间 的 热 漏 ， 系 统 从 
温度 为 Ti 的 低温 热源 吸收 热量 ， 放 出 热量 给 温度 为 
Ty 的 高 温 热 源 ， 如 图 1 所 示 。 该 循环 是 一 个 不 可 道 
量子 斯 特 林 致 冷 循环 。 过 程 1-2 为 等 能 量 过 程 , 在 此 
过 程 中 系统 的 总 能 量 保 持 不 变 ， 由 此 可 得 : 
E, = pubEiit Di 天， 一 
pnbntpwEr2= prbrtpybys (9) 
式 中 Pi,En(i=1,2,L;s=1,2) 分 别 表示 概率 密度 与 对 
应 能 级 的 能 量 本 征 值 ,i=L 表示 任意 态 ，n 为 量子 数 。 
把 式 (4) 和 (5) 代 入 式 (9) 可 得 : 


(4-3p)/B =(4-3p,)/P =(4-3p,)/B 
(10) 


F 4 4(Ea,Ta) 


QH (E313) 


地 ce 
1(E1,T1) 
Qt 2(E2,T2) 
有 


la L2 


图 1 量子 斯 特 林 致 冷 循环 示意 图 . 


Figure 1. Sketch of a quantum Stirling refrigeration 


cycle. 
由 式 (8) 和 (10) 可 以 求 得 在 等 能 量 过 
势 阱 壁 上 的 力 为 : 
F,=Xh (4-3p)/ (mLE) 
等 能 量 过 程 中 系统 的 总 能 量 保 持 不 变 ， 系 统 从 低温 
热源 吸收 的 热量 全 部 用 来 对 外 做 功 ， 因 此 : 


寸 程 1-2 中 ,作用 在 


(11) 


D2.2 
Th 1 daL 
@ = Wi, = (4—3p1) 
12 12 | 而 天 11 L (12) 


=zh" (4-3p mL) In(L, /1) 
同 理 , 在 过 程 3-4 中 ,系统 释放 给 高 温 热 源 的 热 
量 与 外 界 对 系统 所 做 的 功 等 于 : 


zx’h’ 1 
-| 人 -ap 好 
(13) 
=Xh (4-3p)nD) In(L, 7 万) 
过 程 2-3 与 4-1 为 回 热 过 程 ， 理 想 回 热 过 程 中 ， 
在 回 热 器 中 交换 的 热量 分 别 为 : 


Al 六 

0,; = wt 3p11) PE (4-3p3)] 414) 
Azh’ 1 五 

2 = 7[(4—3p11)——(4-3ps)] (5) 
mn 开 了 


从 式 (14) 和 (15) 可 以 看 出 ，Qw3=Q41, 这 说 明 , 理想 
量子 斯 特 林 制 冷 机 可 以 实现 完全 回 热 。 假 设 高 低温 
热源 之 间 存 在 热 漏 ， 热 漏 率 由 下 式 决定 : 


0 =NA ha/(mD) 
式 中 o 为 常数 。 则 每 一 个 循环 的 漏 热量 为 : 


2. = 0 T (17) 
式 中 ?为 循环 周期 。 在 4-1 这 个 回 热 过 程 中 ， 如 果 是 
理想 回 热 过 程 ， 系 统 会 在 状态 1 结束 回 热 过 程 ， 由 
于 内 部 不 可 逆 性 ， 实 际 的 斯 特 林 制 冷 机 存在 回 热 损 
失 , 因为 非 理 想 回 热 , 系统 在 状态 1” 结 束 回 热 过 程 。 
引入 一 个 非 理 想 回 热 因子 yAy>0)， 则 非 理 想 回 热 过 
程 中 的 损失 可 表示 为 : 
Q, = HO (18) 
因此 ， 从 低温 热源 吸收 的 净 热 量 与 释放 给 高 温 
热源 的 净 热 量 分 别 可 以 表示 为 ; 

2 =C.+O+C. (19) 

Q,=0Q,-0,-Q, (20) 

粒子 在 某 一 能 级 上 的 概率 密度 由 吉 布 斯 分 布 

为 : 


(16) 


Pi=1/{1l+exp[z ry AKT.2mE)} QD 


Bi=1/{1+exp[z ry Ax kT, 2mD)]} (22) 
PBi=1/{1+exp[z rn Ax kT 2mD)]} (23) 
Pi=1/{1+exp[z ry AKT2mD)]} (24) 


3， 循 环 周 期 
假设 势 阱 壁 以 速度 "0 运动 ,平均 速度 为 ”可 得 : 
L,—L = | v(1)dt = vr, (25) 


等 能 量 过程 1-2 与 3-4 所 需 的 时 间 为 : 

2=T34=(b -LD)/v (20) 
过 程 2-3 与 4-1 为 等 势 阱 壁 过 程 , 在 此 过 程 中 系 
统 内 能 变化 。 假 设 内 能 随时 间 的 变化 率 为 常数 ， 
为 在 势 阱 宽度 不 变 时 ， 系 统 内 能 是 温度 的 单 值 函 数 ， 
故 有 : 


dT/dt = 7) 


式 中 M 是 一 个 与 系统 能 量 无 关 的 常数 ， 它 只 取决 于 
可 热 材 料 的 回 热 特 性 。 正 负 号 分 别 表示 过 程 2-3(i=1) 
与 过 程 4-1(i=2). 为 简单 起 见 ， 假 设 M/=M2。 因 此 ， 
热 过 程 2-3 与 4-1 所 需 的 时 间 可 分 别 表示 为 : 
ta =(T -TT )/M 
Ti =(T —T)/M (29) 
整个 循环 的 周期 可 以 表示 为 : 
T= 十 723 十 94 十 Zi 二 


2( —L)/v+2(T, —T)/M 
4. 不 可 道 斯 特 林 制 冷 循环 一 些 重要 性 
能 参数 及 分 析 


由 式 (12〉 和 (13) 可 得 循环 的 输入 功 为 : 


! 回 


(28) 


(30) 


xn 1 
W = ln x{(4—3p11) (4—3p31)] (31) 
mL x 


式 中 x=LyL1， 表 示 势 阱 宽度 比 。 由 式 (19) 和 (31) 可 得 
该 不 可 道 量子 斯 特 林 制 冷 机 的 性 能 系数 (COP) 为 : 


cop={(4-3pi)Inx-ar— Hl(4—3p) R={(4-3p,)Inx—a[l2L (x-D/vt+ 


—(4—3ps)/x ]}/{In x(4-3p1) G2) 2(77 —T71)/M]— ul(4—3p1) 
—(4-3p3)/x°]} (4—3ps)/ x}/{2L x-D/y 
联 立 式 (19) 与 (30) 可 得 的 致 冷 率 为 : +2(7 -7.)/M)} 
炳 产 率 为 
1] 元 1 
ee 瑟 (4-3p) 
天 2 
ed | In p;,, 9 
x 2T, x ai | 
1 kmD 1 kmP 1 
十 27 7 (1— Plt A In Pii T2707 Op In Pa 12372 (1— psa)] 
类 给 出 率 为 : 
A 


Ls 
zz DO+G -Dsl= 


[@ DEA{4-3p) nx ar— (4-3p) 
7 mL 


1 
2 (35) 
人 m1 
G3p0/x + 人 2 7(4—3ps)Inxt+ar+ul(4—3p1)—(4— 3p3)/x 1}} 
H 


天 


生态 学 目标 函数 873 为 ; 


A 1 了 人 也 
瓦 = 一 -am =—{(— 和 a {(4—-3p)Inx—-ar—u(4-3p)-(4-3p)/x + -)) 
7 T 全 mL In 


4-3p nxtart (4-3p)-(4-3p,)/0] Ln I 3p11)Inx— 
x tr mb 


1 1 
aor— HI(4—3p11)—(4 eee (30) 


1 Ka 五 
JT 局 (1— p31)] 0 In py +-——s(— pa)]+ 3 a In Pi 二 


( p11) i 
27 ,x a i 27 x 


0 太 允 是 册 下 二 学 目标 必 攻 人 


由 式 (27), (28) 和 (G3) 可 知 ， 对 于 给 定 的 Ts To To Lb 站 15,M = 人 DZX=00L1C 人 =300。 


而 减 小 。 随 着 非 理想 回 热 因子 1 的 增加 ，COP 减 小 。 
图 3 给 出 了 j=0.01， a 取 不 同 的 值 时 ，COP 随 x 的 
变化 关系 。 计 算 中 参数 取 值 与 图 2 相同 。 从 图 可 知 ， 


和 生态 学 目标 函数 都 是 势 阱 宽度 比 x 的 函数 。 图 2 
给 出 了 非 理想 回 热 因子 wk 取 不 同 的 值 时 ， 性 能 系数 随 
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COP 随 着 a 的 增加 而 减 小 ,。 图 4 给 出 了 制冷 率 随 x 的 ” 数 的 取 值 与 图 2 相同 。 从 图 可 知 , 制冷 率 随 x 的 变化 
变化 关系 ， 计 算 中 oa=0.1， J 取 不 同 的 数值 ， 其 他 参 。” 关系 为 类 抛物 线 关 系 ， 由 极 值 条 件 OR/Ox=0 得: 


溉 
[和 D+ -TN 全 -从 IG 3pw) +3p0 rr 7 
2L L Sm 
{4-3p nx-ar- M4-3p) -34-3p)]}=0 
式 (36) 给 出 了 制冷 率 与 势 阱 宽度 比 x 间 优化 关 中 2 te 
系 。 图 5 给 出 了 pc0.01，o 取 不 同 的 数值 时 ， 制 冷 率 上 和 后 
随 x 的 变化 关系 ， 制 冷 率 随 着 热 漏 的 增 大 而 降低 。 | 
图 6 给 出 了 无 量 纲 生态 学 目标 函数 与 制冷 率 之 
间 的 优化 关系 (E*-R)。 从 图 可 知 ，E”-R 优化 关系 曲线 jial 
为 扭 叶 型 ， 存 在 着 一 个 生态 学 目标 函数 的 极 大 值 点 中 
与 一 个 制冷 率 的 极 大 值 点 ， 且 每 一 个 生态 学 目标 函 EE 
数值 ( 除 极 大 值 点 外 ) 对 应 两 个 制冷 率 取 值 ， 显然 要 取 图 3 a 对 COP 与 势 阱 宽度 比 的 关系 的 影响 . 
对 应 的 制冷 率 较 大 的 状态 点 。 热 漏 不 改变 曲线 的 类 


Figure 3. Effects of & on COP with the ratio of 
型 ， 但 随 着 热 漏 的 增加 ， 生 态 学 目标 函数 值 和 制冷 
率 都 变 小 ， 在 生态 学 目标 函数 取 值 较 小 时 ， 热 漏 对 


Potential well width. 


其 影响 较 小 ， 在 生态 学 目标 函数 取 值 较 大 时 ， 热 沁 a tooo 
Ca=0.1.h=| 
对 其 影响 较 大 ,图 7 给 出 了 /对 ER 关系 曲线 的 影响 ， | Seas 


比较 图 6 和 7 可 知 ， 热 漏 对 E*-R 关系 的 影响 比 非 理 | 
想 回 热 因 子 的 影响 大 ， 即 热 漏 是 比 非 理想 回 热 更 重 | 
要 的 不 可 逆 因 素 。 | 
02 
14 
12[ 1.0=0.1,p=0.01 1 15 2 25 3 3.5 4 
2.0=0.1,u=0.1 
1r 3.a=0.1.n=0.3 


图 4. 1 对 致 冷 率 与 势 阱 宽度 比 关 系 的 影响 . 
Figure 4. Effects of hn on cooling rate(R) with the ratio 


of potential well width 


0 5 结论 
1 1 本 文 建立 了 以 无 数 个 无 限 深 势 阱 中 的 粒子 为 工 


质 的 不 可 逆 量 子 斯 特 林 致 冷 循环 模型 。 基 于 薛 定 雇 
方程 的 求解 分 析 了 该 不 可 逆 量 子 斯 特 林 制 冷 机 的 性 
能 ， 给 出 了 性 能 系数 、 致 冷 率 、 烂 产 率 和 生态 学 目 
标 函 数 ， 分 析 了 热 漏 系数 与 不 完全 回 热 因子 对 性 能 
系数 、 致 冷 率 以 及 生态 学 函数 的 影响 。 制 冷 率 与 生 


图 2 W 对 COP 与 势 阱 宽度 比 的 关系 的 影响 . 
Figure 2. Effects of 1 on COP with the ratio of potential 


well width. 
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图 5 a 对 致 冷 率 与 势 阱 宽度 比 的 关系 的 影响 . 
Figure 5. Effects of o on cooling rate versus the 


ratio of potential well width . 
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图 6 a 对 致 冷 率 与 无 量 纲 生态 学 函数 关系 的 影 
响 . 


Figure 6. Effects of o on cooling rate versus 


dimensionless ecological objective function. 
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图 7 4 对 致 冷 率 与 无 量 纲 生态 学 函数 关系 的 影 
响 . 
Figure 7. Effects of 1 on cooling rate versus 


dimensionless ecological objective function. 


型。 


制冷 率 与 生态 学 函数 的 关系 曲线 是 回 原点 的 扭 叶 


存在 着 一 个 生态 学 目标 函数 极 大 值 与 一 个 制冷 


率 的 极 大 值 。 这 两 个 极 大 值 之 间 区 域 是 该 制冷 机 的 


最 优 运 行 区 间 。 本 文 结 


阱 ! 


的 理解 。 


有 助 于 加 深 对 以 无 限 深 势 
粒子 为 工 质 的 不 可 逆 量 子 斯 特 林 制冷 循环 性 能 
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